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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1268 

ISENTROPIC PHASE CHANGES IN DISSOCIATING GASES AND THE 
METHOD OF SOUND DISPERSION FOR THE INVESTIGATION OF 
HOMOGENEOUS GAS REACTIONS WITH VERY HIGH SPEED 
By Gerhard Damkohler 


I. INTRODUCTION AND STATEMENT OF THE PROBLEM 


In seeking to Investigate kinetically those homogeneous gas 
reactions which proceed at very high speed or at very high temperatures 
we are frequently led to the field of isentropic changes of state in 
dissociating gases, a field which as yet has received no detailed 
theoretical treatment. This is the more remarkable since more than 
30 years ago W. Nernst proposed the investigation of very high speed 
reactions kinetically through measurement of the speed of sound at various 
frequencies, that is j a method whereby isentropic or approximately 
isentropic changes of state are impressed on the gas and it is observed 
whether the reactions corresponding to the equilibrium state can or 
cannot follow the sound frequency used, the velocity of the Bound being 
different In the two cases!. This sound dispersion method,_in which more- 
over instead of the sound velocity the sound absorption can likewise be 
measured as a function of the frequency, has subsequently often been 
used, particularly after it was shown by H. 0. Kheser 2 that also the 
transition from the energy of translation into internal molecular 
energy (primarily energy of vibration) can be kinetically investigated 
by sound tests of this kind. However it was still impossible to apply 
the sound dispersion method to cases of more complicated reaction kinetics 
since the evaluation formulas of such sound tests always assume a 
knowledge of the kinetic laws3 (in particular the dependence of the 


* n Isentropische ZustandB&iderungen in Dissoziierenden Gasen und die 
Methods der Schalldispersion zur Untersuchung sehr Schneller Homogener 
Gasreaktionen." Zeitschrift fftr Elektrochemie, Bd. 48, Nr. 2, 1942, 
pp. 62-82. (The conclusion of this paper, sections V to VIII, has been 
published as NACA TM 1269.) 

^See dissertation by F. Keutel conducted in the Nernst Laboratory, 
Berlin 1910. 

2 H. 0. Kheser, Ann. Physik 11, 761 (1931) . • 

^See A. Einstein, Sitzungsber . Preuss. Akad. Wiss. Berlin, 1920, 
380 to 385. E. 0. Kheser, Ann. Physik 11, j6l (1931), and H. 0. Kheser 
and 0. Gauler, Physikal. Ztschr. 37 , &77 to 684 (1936). 1 
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reaction velocities on the concentrations) and only the remaining 
undetermined reaction velocity constants are obtained. Thus the problem 
still remained of deriving an evaluation formula for such sound tests 
which is free from any restricting assumptions as regards the unknown 
reaction velocities. 


The knowledge of the isentropics .of dissociating gases is also 
required, however, for the understanding of numerous gas dynamic 
phenomena in hot combustion gases. Such phenomena can be very well 
studied in pipes for which purpose several relations of G. Damki'hler and 
Adam Schmidt were recently obtained. ^ In these relations there again 
always occurs the isentropic exponent which for nondissociating gases 
agrees with the ratio of the heat capacities Cp/Cy = k. This, however. 


does not hold for dissociating gases. But in this case too it was 
recently found by G. Damki'hler and R. Edse for several examples in 
as yet unpublished computations that the isentropic can be represented 
by the general polytropic equation pV m = const, (p = pressure, Y = volume), 
even over large pressure intervals. If, for example, the gases of the 


stoichiometric CO combustion (CO + i O 2 as the initial mixture) which 

are assumed initially to be at 40 atm. and 3450° K are expanded 
isentropically to 6 atm. and then again to 1 atm.5 ? the temperatures 
assume the values 2973 and 2615 ° K; the total number of moles changes 

from 1.197 to 1.142 to 1.099 (referred to 1 mole CO + ^ mol 0 2 as 

initial mixture) so that the isentropic exponent m becomes 1.115 in 
the upper pressure interval and 1.103 in the lower pressure interval. 

The value of m thus changes by only about 1 percent even in this large 
interval . 


In the following we can therefore restrict ourselves to the 
treatment of the differential isentropic changes in state in 
dissociating gases. It will be assumed, however, that several 
dissociation equilibrium states occur since this is practically always 
the case for high-temperature reactions. In the dissociation of pure 
COg, for example, we have to take into account both 

2 C0 2 l=±2 CO + 0 2 (1) 

4;Damk6 , hler, G., and Schmidt, A., Ztschr. Elektrochem. 47, 547 to 567 
(1941). 

5 Such pressure relations occur for example in the case of detonation 
heads in pipes . 




NACA TM 1268 


3 


and 


0 2 <=»2 0 


( 2 ) 


In the combustion of hydrocarbons with oxygen we have to consider 
three additional equilibria 


2 Hg0?=&2 H 2 + 0 2 , 

(3) 

2 E 2 0T ± -‘S 2 + 2 OH, 

W 

H 2 «±2 H 

(5) 


and in the presence of nitrogen, also the reaction 


N 2 + 0 2 -?±2 NO 


(6) 


In what follows the problem of the differential isentroplc changes 
of state in dissociating gases will be treated quite generally, namely, 
the attainment of complete equilibrium at each instant on the one hand 
and on the other the incomplete equilibrium state in the sonic field. 

A generally valid formula for the reaction-kinetic evaluation of sound 
dispersion tests is obtained. Almost all physical effects which like- 
wise give rise to sound dispersion and can therefore act as disturbing 
factors are discussed and quantitatively estimated. General points 
of view are obtained for the experimental carrying out of sound tests 
in the high- temperature range. 


II. THE DIFFERICNTIAL ISENTROPIC EXPONENT m FOR THE CASE OF 
COMPLETE ATTAINMENT OF CHEMICAL EQUILIBRIUM 


We consider a definite quantity of an ideal gas mixture of 
various kinds of molecules J, among which exist n thermodynamically 
mutually In dependent equilibrium conditions. To these correspond n 
overall reactions, which however can be arbitrarily chosen to some 
extent. 
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T 

V 


The following symbols will "be used: 

absolute temperature (°K) 
volume of gas considered 

number of g-moles of the type j contained 
in Y 



M J 


total number of g-moles contained in Y 
molar weight of the particles of type j 



density 


°j - V T 
p J 3 

p 

E 


g-moles of the type J in unit volume 
partial pressure of particles of type j 
total pressure 

gas constant = 1.&86 cal/mole deg = 8.313 x ID 7 
erg /mole deg 


-PJ 


molar heat at constant pressure of the 
molecules of type j 


£vj ~ £pj " E molar heat at constant volume of molecules 

of type j 

H enthalpy of gas system considered 



partial molar enthalpy of molecule of type j 



heat capacity of the gas system considered 
at constant pressure and constant number 
of moles Nj 

heat capacity of the gas at constant volume 
and constant number of moles 


m 


exponent in the isentropic pV™ = constant. 
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The 1th overall reaction is assumed, to satisfy the stoichiometric 
equation of the type 


V 11 X 1 + 


52 


*5 + | v 6 i ]*6 


. + W 7 


(7) 


■where 


X J 


is the chemical symbol for the molecule of 
the type 3 



the absolute value of the stoichiometric 
conversion coefficient of the molecules 
of type J (first index) for the 1th 
reaction (second index); for the reaction 
going from left to right, corresponding 
to the Eucken notation^ we have 


V J2 

V J2 


> 0 

if 

the molecule 

of type 

J 

vanishes 

< 0 

if 

the molecule 

of type 

i 

originates 

= 0 

if 

the molecule 

of. type 

‘3 

remains unchanged 



*1 = 


heat of reaction at constant pressure for 
the 1th chemical reaction 


A differential isentropic change in state, indicated by 5, in 
the gas mixture considered must on the one hand satisfy the first 
fundamental law, i.e. 


V5p = SH 



5T + 

P N 3 



8p + 



°Precisely the opposite sign convention is used by American 
investigators. 


( 8 ) 
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or with the relation 






'/Tffj 


with the above notation. 


= 0 which holds for ideal gases and 


J- 


CpST - VBp + 2_ H J 51T J = 0 


(9) 


On the other hand each of the n independent overall reactions must 
satisfy the law of mass action, which for the Zth reaction assumes 
the form 


j 

^~V jZ lnpj = lnK z (10) 


where K 1 denotes the equilibrium .constant which depends only on the 
temperature. For the differential change in state we have therefore 
n further determining equations: 



( 11 ) 


In equations ( 9 ) and (ll) we now replace the differential changes 
8 T, 8 pj, and BNj by the magnitudes SV, Bp, and 5z m where 5z m 

in dicates how often the mth reaction equation ( 7 ) for the differential 
change In state has proceeded from left to right. We have 

m 

<®J - (12) 


BN 



v m^ z m 


j 


m 


I 


W-mBZm 
n m ^m 


( 13 ) 

(14) 


From 
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follows 


ST Bo , BY 5U Sp SY , 5z m 

T "p Y“N“p Y Z_ Vm HT 


(16) 


and from 





(17) 


there is obtained. 


m 

+ Y" ( _ A^m 

pj p h Hj - p Z— r m J m nJ ht 


(18) 


Substituting these egressions in the equations (9) and (ll) and using 
the relation 



there is obtained as the expression for the first fundamental law: 


( 20 ) 


and the law of mass action of the Zth reaction is given by 

I? *(5 - A *t (A *t r 3 - - ° <-> 
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For briefness we write: 

* 

A-l ~ KV l ~ [ K “ 'O RT (.22.) 

W ^ 

Blm = — v m + 2_V J1 V |}I11 - V 2 V m (23) 

3 

c l = ET " v l ■- ( 2j+ ) 


hence 


w l 

HE = A * + * C 2 


(25) 


and obtain from the 1 + n equations (20) and (21) the system of 
equations 


p 5V 


Al 


C 2 


p 5Y 


C n ^ m + B r 
n p 6Y * 


8 z-l 

Y 

it 

8Y 

5 zq _ 

Y 

- IT 

5Y 

Sz^ 

Y 

■ IT 

8Y 

5z-j_ 

Y 

■ N 

5Y 


5Zp y 

A_ —l— 4 . 

■^O TVT CTT ' • 


bz 2 J_ 

- IT 5Y 


Szp y 
®22— 5Y 


5z 2 Y 

+ B„ 0 -rr- rr: + 


, A &Z n V 
' + 1 57 ‘ K 


+ B 


5z n t _ 


N 5Y 


= "(Al + KCq) 


+ ^ = “( a 2 + kC 2) 


>(26) 


Sz 


n Y 


• + Bnn jj gy - "(An + K C n ) 
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From the above system of equations the isentropic exponent m can he 
imme diately computed. For, from 


pV 11 = const . 


( 27 ) 


we have 



&£ 

P 

5Y 

+ m ~ 

= 0 


therefore 

0 

A i 

Aq ■ • 

■ ^ 


A 1 

hi 

b 12 ‘ • 

• B ln 


A 2 

• 

b 21 

®22 • • 

« 

• ®2n 
• 

-5p Y 

m = — * ~= = k + 
p oY 

• 

An 

■ 

• 

®nl 

» 

Bn2 - ' 

• hm 

1 

A 1 

A2 • ■ 

• An 


Cl 

B 11 

b 12 • ■ 

• B ln 


c 2 

B 21 

Bg2 • • 

• ®2n 


Cn 

Bnl 

Bji2 • • 

• Bnn 


(28) 


= K + 


& 

*2 


(29) 



The isentropic exponent m in dissociating gases is thus referred, to the ratio K of the 
molar heats and the ratio of two determinants Aj_ and Ag. The dimensionless magnitudes in 

the latter Aj, B^, Cj are immediately given as soon as the composition of the gas and the 

heats of reaction of the n overall reactions are known. In order to avoid any errors in 
writing the indices it is test to prepare a table of the n overall reactions as we shall see 
below. The two dete rminan ts Ap and &£ are to a large extent symmetrically constructed and 

contain 1 + n rows corresponding to the first fundamental law and the n mass action laws for 
the n thermodynami cally independent chemical reaction possibilities. It is immaterial in what 
manner the latter are expressed. If, for example, the original overall reaction (2) is replaced 
by a reaction (?) which is obtained by the addition of the reaction (l) taken x Mmes and 
the reaction (gj taken y times there follows from (22) to (24) 

Ag = k(xVi + yV 2 ) - [k - l] ~ X ^ 1 R ^ - yW2 = xAp + jA 2 (30) 

xW-r + yWo 

^2 - pm (*Vp + yVg) = xCp + yC 2 (31) 


■^ur_ 4 . tAj - 


Ma +y (xvj! + yVj 2 )Vjm jjj -UVl + yv 2 )\k = xBiju + yB^ for (g) / 


m 


(32) 


u 

*12 - + J’M?) + / - n (xy,. + JVT.) - **u + 7*12 for © J @ 

J (33) 


< 


> » 


*■ r 


h 
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Ti_r _ 


3^-1 + 

i 

RT . 


± tr*» i_ \ 

1 j vy>i 


• y ** y jj- 


^2 b _ 

J V J2 / H, ^ 


»l_ _l_ tt**i ^ — ■v2lJ_ __ J_ 'wffTL . _L trwTl^ > 

' j y 2 / - ■*■ ^ i i 1 * n v j - , j _2 1 


-rr 2 n . 
j 


(3^) 


Ve then ha,7e fcsr the detsiTnlnstion In the denossluator in ( 29 ) 



i 

A 1 


’ ‘ \ 


1 

A 1 

(xA x + jA^) 

... An 


Cl 

B n 

b 12 ■ 

• ‘ ®La 


C 1 

B n 

(*%! + yB^) 

■ * • ^ 

= 

C 2 

■ 

# 

• 

• 

^22 • 

■ 

• 

♦ • B 2n 

• 

• 

1= 

(xCtl + yC 2 )(xB 11 

■ 

+ JZqJ 
• 

• 

( x2b h + V 

1 + y^2i + y 2B 22 1 

\ f 

... ( xB-) ^ + yBpn ) 
• 

■ 


• 

n 

'■'n 

n 

-°nl 

Tk — 

^n2 • 

• 

• ■ -Rnn 


• 

r* 

• 

ti 

^nl 

/ Tk ■ Tk \ 

\- UD nl T y- D n2-' 

Tk 

. . . -Pnn 


(35) 


P 
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last d e t e rminflnt remains unchanged, if W e subtract the second, row 
multiplied, by x from the third, row. There is then obtained: 



1 

Al 

(xAi + yAg) . . 

• An 


C 1 

B n 

(xBn + yB 12 ) . . 

• B ln 

Ag = 

yc 2 

ys 2 i 

y(xB 21 + yB 22 ) . . 

• y B 2n, 


C n 

2^1 

(xBni + yB^) . . 

* ®nn 


If we now subtract the second column multiplied bv 
column there is obtained: 


x from the third 



1 

A 1 

y&2 

• • ’ An 


c i 

Bn 

y £ i2 

• • • -^ln 

^2 = 

yc 2 

• 

7 B 21 

CVJ 

• * • y ®2 n 



B nl 

*^2 

' ■ B nn 


(37) 


Taking the common factors 
the determinant we ha ve 


of the third row and the third column outside 


^2 = (38) 

In a similar manner we have for the new determinant in the numerator 

4 - y 2 Ai (39) 
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so that the quotient "becomes 


^ 

A 2 ^2 


( 40 ) 


Thus it has been mathematically proven that the expression (29) 
for the isentropic exponent m is independent of the manner in -which 
the n thermodynamically independent reaction possibilities are 
expressed by n overall reactions. Physically of course no other 
result could have been expected. 


III. THE INCOMPLETE ATTAINMENT OF CHEMICAL EQCILIHRItM IN THE 
PLANE SOUND WAVE (NEGLECTING THE EFFECTS OF HEAT CONDUCTION, 
FRICTION AND DIFFUSION ON THE SOUND DISPERSION) 


If a sound wave is allowed to pass through a gas the latter is 
subjected to periodic compressions and expansions which, to a first 
approximation, are isentropic. In the case of a dissociating gas for 
sufficiently high sound frequencies the chemical reaction velocities are 
no longer sufficient for the' attainment of complete equilibrium and we 
have the case of sound dispersion. In the theoretical treatment of the 
latter, which is partly also a problem of flow theory, three conservation 
laws are to be taken into account: the conservation of masB (continuity 

equations), the conservation of energy (first fundamental law of 
thermodynamics) and the conservation of momentum (Euler equation of 
hydrodynamics) . 

We consider a definite gas mass, i. e. a volume V, fixed with 
respect to the flow, of such small dimensions that the same state may 
be assumed throughout. This volume changes with time by differential 
amounts in the following manner: 

Continuity equations: 


Hi. 

dt 



v j: 


dz m 
m dt 


(41) 
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dN 

dt 


m 


d£m 

m dt 


1 dv; 

V dt 


6 w 


m 

(43) 


where w is the flow velocity at a point fixed in space over which a 
plane sound wave passes in the positive or negative direction. 

Energy equation: 


dE _ dJ 
dt “ " p dt 


(44) 


or on substituting the enthalpy H = E + pV and after some computation 


dT - 
dt 


V 

dt 





(45) 


Momentum equation: 



(46) 


Equations (4l), (42), and (45) agree with the previous equations 
(12), (13), and ( 9 ). Equations (43) and (46) are respectively the hydro- 
dynamic continuity and momentum equation. In equation (46) the added 
friction term is neglected and in equations (44) and (45) the added heat 
conduction. The particle interchange through diffusion between the 
neighboring gas elements is likewise not taken into account. 

For a sufficiently small sound intensity all the state variables 
carry out sinusoidal vibrations about corresponding mean values 
characterized by the index zero. To a first approximation they agree 
with the corresponding magnitudes of the equilibrium equations of the 
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medium without sound waves, 
p = p Q + Bp 
V = T q + BY 
T = T 0 + ST 
p = p Q + Bp 

H j - N Jo + W J 

Iff = IT + BIT 
0 

C J = C J 0 + 5c J 


7 We may therefore set: 


Sp = Re ( 

^p e gs + ktj 

BY = Re | 

[vegx + 

BT = Re 1 

j^Tegx + ktj 

Bp = Re ( 

' re gx + ht ) 


r >■ 


SNj = Re ^Nje 6 * + h ^j 


8H = Re IffeS* ++llt 


*M = Ee £j e^ + « 


5z m = Re ^ z m eS^ + ktj 
Sw = Re fwe® 1 + 


( 47 ) 


■where Re() denotes the real part of (). The underlined symbols are 
complex magnitudes which take into account the phase angles between 
the various magnitudes and are assumed independent of x and t; g is 
complex and h is pure imaginary. They depend on the frequency f, 
the sound velocity a and the damping constant 7 for the sound 
amplitude in the following manner: 


gx + ht = 2itf i t + $ T- 7 x 

s 


m 


^The difference between the mean pressure value pq in the sound 

wave and the pressure in the medium without sound is, for example, of the 
order of magnitude of the sound radiation pressure. The latter is 
proportional to the energy density in the sound wave and therefore 
proportional to the square of the sound amplitude. The above mentioned 
pressure difference can therefore, for sufficiently small sound amplitudes, 
be neglected as a small second order magnitude. 
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or 


a 


2rtf 



and 



= imaginary part of g, i.e. Tm(g) 
a — 

i2?rf + 
a 


^ and y - real part of g. 


since gx = -f 



( 49 ) 


(50) 


For the time variation in moving forward with the flow along the path 
curve dF/dt of any magnitude F in the equations (4l) to (46) we 
can also write 


dF _ 8 f + w dF 
dt 8t 8x 


(51) 


where dF/St and 3 f/8x are the partial derivatives at a point of the 
spatially fixed coordinate system, relative to which the gas particles 
move periodically with the velocity w. In equation (51) the second 
term on the right vanishes with decreasing sound intensity more rapidly 
than the first term on the right because the former is proportional to 
the second power while the latter is proportional to the first power of 
the sound amplitude. We may therefore throughout replace the 
operation d/dt by d/dt and use it directly on the expression (47). 

The magnitude dz m /dt in the initial equations (4l), (42), and (45) 
indicates how often the mth overall reaction occurs in the volume V 
of the gas and Is thus proportional to the mth reaction velocity. It 
depends naturally also on the velocities of the individual elementary 
reactions hut these do not enter separately in the final formulas to he 
derived for the sound dispersion, since they are not contained in' the 
initial equations (4l) to (46). Hence only the velocities of the over- 
all reactions can primarily be determined from sound dispersion 
experiments and not the velocities of the elementary reactions. 


Added hy the MCA reviewer for clarity. 
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We consider again the Zth total reaction 

| v iz| x l + I v 2z| *2 + * * • f^| v 5Z x 5 + v 6 z | x 6 * +•'••+ ( 7 ) 

and denote by 

t?2, its frequency from left to right per unit of time and volume 
at a given instant j 

tf 7 its frequency from right to left per unit of time and volume 
at a given instant, 

Uj its corresponding frequency from right to left and left to right 
in the equilibrium state ( state in the medium without sound = mean value 
of the state in the medium with sound), where = tfj = ^z • 

For dz^/dt there is then obtained 

M,V( Uj -Uj) (52) 

at 


1 

The reaction velocities and depend on the temperature T 

and the various concentrations cj. Since only very small changes in 

jBtate arenas sumed impressed on the gas by the sound wave the magnitudes 
and Uj differ only by a small amount from the reaction "velocity Uj 

in the sound- free equilibrium state. We may therefore develop the first 
order to arrive later at more general formulas free from special 
assumptions: We start not with but with In tfj: 


In Uj = In Uj 


+ l- SA 8T + Sc 

3t - 8c . j’ 


— ) ^ —i J A 

U 7 8 In U 7 v — 8 In U 7 

In _L = L 5T + > L - 5c, 

U z 3T 3c 1 


Uz - Uz\ U! - \Ji 


= In 1 + 
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and similarly for the right to left reaction 

i 


(r~ 

n u 2 

In — = 
U 1 


3 In Uj 

St 


5T + 


S In U2 

Sc 4 


Sc j ~ 


In 


(r- 

Uz 


Uz\ 


U2 


U2 


- Uz 
U2 


(5^) 


In the neighborhood of the equilibrium state the velocities U2 

and U7 are no longer entirely independent of each other but mathematically 
L **» 
must possess such form that in the state of equilibrium itself where 

and become equal the mass action law applies, i. e. for our 

reaction (7): 



const. 


(55) 


This is possible only If we set: 


U2 = kzci 1 


*11 


C2 1 


F(T,ci,C2, . . .c j) 


(56) 


and 


e- 



*z c 5 




F(T,ci,C 2, • • -Cj) 


(57) 


In both equations on the right there is the same factor F(T,ci,C 2 , ■ . .cj) 

which otherwise may quite arbitrarily depend on the temperature T and the 
concentrations Cj. It is not difficult to see that by means of the 

expressions (5 6) and (57) all reaction-kinetic laws for the reactions in . 
either direction can be represented, subject to the only restricting 
condition, namely, that in equilibrium they must lead to equation (55)* 
The magnitudes I?2 and ^2 in equations (56) and (57) are the reaction 

velocity constants which depend only on the temperature corresponding 
to the Arrhenius equation 
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and 


5 In kj 

8t 


(- 



(59) 


where q^ and q^ are the heats of activation for the reactions in the 

corresponding directions. Substituting the expression (56) to (59) in 
the equations (53) and (5*0 and then subtracting the two latter from 
each other there is obtained after brief computation 


— ) 
U, 


- u z 


-4 (- 

= * 1 : ^ 

RT 


5 In T 




v 3l 5 


In c . 


(60) 


The as yet undetermined factor F(T,c-pC£, . . .c^) in the expressions (56) 

and (57) dropped out in the subtraction. 8 The stoichiometric coefficients 
in equation (60) again have the signs corresponding to the Eucken 
notation given above. In equation (6 0) there occurs the difference 

in the heats of activation q j - q^ which may be expressed in terms 


of the heat of reaction . 


^c 10 

v iz 

n 

ro 

0 

l V 

l 2 1 

• ■ • 

V50 

v 52 

c 6o 

v 6l\ 

• 

• • • 


From the kinetic equilibrium condition 

- y | ^ 

= 1, or. in k* - In k z = - > v^c l0 (6l) 


and the thermodynamic equilibrium condition 

J 0O C J0 + in R + In t] 


In Kj = ^_Vjjln p 


(62) 


we have 


— > 4 — 

In kj - In kj = 


- In K z + Vjln R + Vjln T (63) 


^At this point it is readily seen why we developed above the Taylor 
series not for Uj but for In Uj 
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or differentiating with respect to the temperature T with account 
taken of ( 58 ), (59), and (ll) 


-» <- 
9l - = 

RT 


ll 

RT 




(64) 


Finally making use of 


Slncj=-SlnV+81nRj 


-SV 

V 




_^m 

N 


( 65 ) 


and the previous relation 


8 In T 


, 51 + 5p 'y- v ffm 

V p ki_ m N 


( 16 ) 


there is obtained from equation ( 60 ) after brief computation 


Uj - Ul 
Uz 




( 66 ) 


The above equations (66), of which there exist as many as there are 
thermodynamically independent overall reactions, completely correspond 
to the n previous equations (21) as is seen by comparison of the 
factors in parentheses. Equation (66) now contains also, however, the 
reaction velocities Uz, tfj, and Uj of the 2th total reaction, as 
must be the case in the sound dispersion field. For very small 
frequencies f or very high reaction velocities Uj equations (66) go 

over into the previous equations (2l) . For, from equation (52) there 
is obtained, using equations (47) and (48) ■ 


*4 

Uz 


Uz 


u. 


_ _!_ £51 - 
U z v dt 


Re< 


2rtfiN — Z = gx + htl 


u z v 


N 


( 67 ) 


an expression which is proportional to the sound frequency f and 
inversely proportional to the reaction velocity in equilibrium. 



I 


If we Introduce the dimensionless magnitudes 


q> 77 = im = l nf ^ c x 
11 U Z V U 2 


and again use the previous notation 


ftyu i P/niflt.innR f 1l 7^ ■f.liArA la n"h+.ti1 ■Prmn omio+. •! rvrx ( )l.R ^ trVii 0 V 1 Iqd/Io +a +Vio nne 

— -i \ ■ j / «**»* «*- *-> u a ■». wAu, aua* \ ^ r f y nij.xvu -*-*-» uu,u ww vuw j/A >w t a wiau 

equation (20), and from equations (66), after brief computation, the system of equations 

A ^ 

"n h 


v p 

K~ + = + 

V p 


a ^ 

Al ¥ 
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3o 


V p 

( A 1 + K C Vy + C X ~ + (B lx + i(Pn)^ + B 12 Y + 
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In 


If 


= 0 


> 


V P Z-, Zp 

(Ag + ItCg)— + Cg — + Bg-L — + (Bgg + ic P22)~ + 

■ V . p . N JN 


In 


N 


Bpn — = 0 

> N 


( + KC n)~ + c n _ + Bnq 

V p - 1 - 


zi Zg 

~ + Bno — + 

^ If 


+ (Bnn + Kp^n) ~ = 0 


If 


(68) 


. .. r -Tl W 2 

^ =KV l - L K - ^KC 

Vtt) 

B Zm = +y~ v .i Z v .1m ~~ ~ v Z v m 

( 23 ) 

W 7 

n-, = — s. - Vt 

t KT <- 

fpu 

\ — ' f 


( 69 ) 
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From the above there is obtained for the ratio of the relative pressure cha n ges to the relative 
density and volume changes respectively, in the sound wave 



U 

A 1 


A2 • . - 
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A 1 

( B ll + i<Pll) 
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B ln 
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^1 
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pr pv K + . 
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(70) 


ro 
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Equation (70) goes over into the previous equation (29) if* the phase angle^ in the 
expression on the right vanishes, i. e. if in the determinants Ai‘ and £ q ' the 


Qriv. 

~ xuB puetse tmgJ.e 








y iOiaxccibtits uy nuw muon biic px<rb blll k outiugts op 


the density and volume changes 8p and 5V respectively- 
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imaginary terms icpq^ vanish. According to equation (68) this is 
the case for sufficiently small sound frequencies f or sufficiently 
high reaction velocities Uj. Thus we have also derived our previous 
equation (29) "by a purely kinetic method. 

The determinants Aq‘ and A q' can he obtained from the 
determinants Aq and Aj respectively of equation (29) by substituting 
in the latter for the Bi m values those values which satisfy the equation 


B Zm' = B lm + i( ?Zm 


(71) 


with the added conditions 


Vim 


fgtrfN 



for l = m 
for l / m 


(72) 


The ratio pp/pr given by the expression (70) can also be obtained 
from the purely hydrodynamic equations (43) and (46) by substituting 
the expressions (47) and replacing, as before, d/dt by d/dt. We then 
have from (43) 



and from (46) 


phw = - gp 


(73) 


(74) 


and multiplying the two last equations 


h 2 r = g2p 


or finally 


S = 


Th 



(75) 


Substituting in the above equation the expression (70) and setting 

p = Si = IE 

P NM 


( 76 ) 
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where M is the mean molar weight there is obtained 



Thus the sound velocity a and the damping constant 7 are referred to 
the magnitudes m and cp computable from equation (70). For we have 



The damping of the sound wave per wave length l = a/f is given by 

yl = 2rttan & Jg ircp^l + (80) 

Since physically the damping constant 7 enters only as a positive 
value in the expression ( 48 ) it follows from (7 9) that also the phase 
angle cp, which always amounts to only a few degrees (see below), is 
positive. Hence according to equation ( 70 ) the pressure change 8p 
leads the volume and density changes SV and 5 p respectively. It has 
already been shown clearly by H. 0 . Kneser (Ztschr. Techn. Physik l6, 
pp. 213-219 ( 1935 )) that only in this way is there obtained for periodic 
changes of state (on the pV-diagram) a closed curve leading to energy 
absorption. 

Equations (70) and (78) and ( 79 ) can be ulitized in the following 
manner for describing the kinetic relations in multiplying dissociating 
gases: For a constant gas composition, fixed pressure and temperature, 

definite numerical values are assumed for the n reaction velocities Uj, 
and m and cp values are computed by equation (70) for various sound 
frequencies f and from these with the aid of (78) and (79) the sound 
velocity a or the damping constant 7. There is then determined, in 
a series of measurements for the same gas composition, the same pressure 
and the same temperature, the sound velocity or the sound absorption 
as a function of the frequency. Finally the experimental a-f and 7-f 
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curves are compared with the theoretical curves and it is seen which 
expression for the assumed U7 values gives best agreement with the 
test curves. If, in particular, of the n reaction velocities the 
numerical values of n - 1 are known (from other considerations), the 
last unknown reaction velocity can be given exactly for the chosen 

test conditions. By systematic variation of the latter in various series 
of tests in which again the total frequency range is traversed the last 
reaction velocity U n can be determined as a function of the concen- 
tration Cj and the temperature T, provided that these relations are 

already known for the remaining n - 1 reactions. 10 Hence no special 
reaction kinetic law is a priori assumed for U n but the latter can 
be quite arbitrary and is determined by comparison of the experimental 
and theoretical curves. The new procedure described for the evaluation 
of sound dispersion tests thus differs fundamentally from the procedure 
used thus far since the latter make special assumptions on the kinetic 
laws of all the individual reactions whereas we now require only a 
knowledge of the velocities of n - 1 total reactions. This will be 
seen more clearly be^-ow. 


IV. EXAMPLES FOR THE APPLICATION OF THE NEW FORMULAS 
1. Older Cases of the Reaction-Kinetic Sound Dispersion 


Before we apply the newly derived equations on the problems of 
particular importance of high-temperature 'dissociation it will be shown 
that naturally all the cases of sound dispersion already considered by 
other authors are correctly described by the new formulas. • We shall 
restrict ourselves to the case of Incomplete attainment of an individual 
dissociation equilibrium, which was first computed by A. Einstein, 
and furthermore the incomplete attainment of an energy term (vibration 
term) in a molecule, as investigated by H. 0. Kneser.12 

Einstein considers only a single dissociation equilibrium in which 
from a molecule Xq there arise two equal molecules X2 so that the 
reaction equation reads: 


Xq £=> 2 X 2 + W x (81) 

We shall denote it by the ’reaction index' 1. 

^This is also superfluous if the sound dispersion is "mainly" 
brought about by a single total reaction. See the example of the C0 2 
dissociation below. 

-^Einstein, A., Sitzungsber. Preuss. Akad. Wiss., Berlin 1920, 380 to 385. 
■^Kneser, H. 0., Ann. Plpysik 11, j6l (1931) 
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For the doubly indexed stoichiometric coefficients^ there is obtained 


'll - 1 'fel = - 2 


and therefore 


Vl = 



V jl = “ 1 


(82) 


Einstein does not make use, as we do, of the heat of reaction W]_ at 
constant pressure but the dissociation heat D at constant volume so 
that ve obtain 


Wi = - (D + RT) 


(83) 


For the constants (22) to (2k) there is therefore obtained 


Al = KM - ( K - 1)^ = - K + ( K - 1 )/^ + 1 ] = ( * - l)j| - 1 (8k) 


J 


Bll = Ulyj, + vu 2 JL . Vl 2 = (-2. + 1 + (JL + M . 1 = JSL + 2L + Jffi 
RT J Nj 1 ^RT j N 2 j RT % N 2 J 

(85) 


W, 

Ci = — - Vi = 

x RT 1 


-2- + 1) + 1 = -=2 


RT 


RT 


( 86 ) 


For the dissociation Einstein assumes a reaction of the first order 

and for the recombination one of second order. This gives in equation (68) 


— > *— p 

Ui = kici = kqc2 


■'■^The type of molecule is denoted by the first index and the 
reaction by the second. 



< 


or 


<J>H = with. o> = 2irf 

*l»l 


(87) 


From equation (70) there is then obtained 
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£ = £ me i9 = P , 
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► = P J 

1 ! .(« - i)[ B n + iq>n - A i c i] - A i 2 

r' p p 

Bn + i<Pn - A]Ci 
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B 11 + i<Pll - A1C1 


(88) 


Since (it - l) C-, + An = 


1, we have r 

(1* - 1) 


E p 
1 “ P < 


1 + 


N 4n\ 


iaiN 


w + \*i + 15 




(« - 1) 
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v_ 

RT + \s 


B- + jfflU 

^2; 


ioatf + JD_ 

t A ™ 


l (K - - 11 


(89) 


or 
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Taking into account the previously given relation R/(k - l) = C v /N 
there is obtained finally 


g 

r 



(91) 


This equation is entirely identical with the expression (19) given in 
the Einstein paper (loc. cit.) if, we note that the term in parenthesis 
(l + 4N]_/N2) is there erroneously written as (l - 4N]_/N2) • 

It may be remarked that the expression (9l) even for the simple 
dissociation equation (8l) is already quite complicated. In particular 
it is difficult to see how it can be generalized when there are several 
dissociation equilibria under consideration. By way of contrast the 
new expressions (29) and (70) of the present paper are of very symmetric 
form and such that the physical meaning is clear. 

M. 0. Kheser 1 ^ considers a gas whose molecules possess a single 
internal energy term and which in the sound field is incompletely excited. 
In order to be able to treat this case with the formulas of the present 
paper we consider the gas as built up of two types of molecules 
and X 2 : the former is assumed not to be excited in the energy term 

under consideration while the latter is assumed to possess the exciting 
energy e per mol. The excitation may then be considered as a chemical 
reaction^. 


X x >X 2 


(92) 


which occurs without change in the total molecule number with the heat 
of reaction 


Wl = - 6 


( 93 ) 


^^Kneser, M. 0., Ann. Physik II, j6l (l93l) 

^As a single reaction it is denoted with the reaction index 1. 
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The stoichiometric coefficients are then 


J 

= l v 2 ]_ = - 1 Vj_ = v ji = 0 (9 1 *-) 

In the computation of k only that part C a of the molar heat C v may 
naturally be used which, independently of the sound frequency, is always 
fully excited (as for example translational and rotational energy) and 
corresponds to the so-called external degrees of freedom: (.’• the index a) 
Hence 


C a + R 

2a 


(95) 


The ratio of the internally excited to the unexcited particles is given 
by the Boltzmann law 


52 

Nl 


-g 

RT 

e — and 



-g 


1 + 


RT 
e — 


thus 


N 

n 2 


1 + e RT 


(96) 


for complete equilibrium the part of the molar heat C v due to the 

vibrational internal energy term is 



2 


( 97 ) 
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after "brief computation for our expression ( 68 ) a"bove 


cPn = 2itfp 



- 2 nf 0 



(103) 


Substituting the expressions (99) to (101) and (103) in equation (70) there 
is obtained 


| = f 

r p 


£ < 

P 


R(l + i 2rtf0) 

1 + C a ( 1 + i 2«f3) +Ci 


or using equations ( 78 ) to ( 80 ) 


1 /mRT 

a ~VT 


7l ~ TtCp 


^RT / R 

C T + C a (2TCfp) 2 ]\ 

5\ . 

lv 2 + Ca 2 (2rtfp) 2 y 

rtfcp 


a 



R £2 2n 2 f2(3 

a[c v (C v + R) + ^(Ca + R)( 
) 


RCi 2rt 2 ff3 


C V (C V + R) + C a (C a + R)(2rtfp) 2 


(104a) ) 


(104b) 


(104c) 


(l04d) 


These equations are entirely identical with those obtained by Kneser. 
Equation- (l04d) represents as a function of the frequency f a bell- 
shaped curve whose maximum lies at 


_1_, f c v (C Y + R) 

2*3VCa(£a + R) 


(I04e) 
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and whose half value width in the case C a 8 C v lies at about 3*8 octaves. 

The example considered above of the incomplete molecular excitation 
in the sound field naturally represents a very simple case. Considerably 
more complicated cases have already been both experimentally and 
theoretically investigated. 17 Of interest is the more or less mutually 
independent excitation of the various natural vibrations in the same 
molecule. A theoretical contribution to this problem has recently been 
given by Kl. Schafer. 18 It is not impossible that in this case, too, 
the formulas derived in the present paper can be successfully applied 
if it is .desired to make no assumptions of any kind on the kinetic laws 
of the vibration excitation and yet derive definite numerical values 
for the reaction velocity. If, for example, each gas molecule possesses 
two different natural vibrations then, in a correspondingly chosen 
temperature range, the molecules could be divided into the following 
types: 


1. molecules Xj_ without any vibration excitation 

2. molecules Xq in which only the first natural vibration is excited 

3. molecules X3 In which only the second natural. vibration is excited. 

Among these three types of molecules there exist two equilibrium 
conditions so that in equation (70) we have 3-rowed determinants. This 
example will, however, not be further considered here since this would 
exceed the limits of the present paper. We shall discuss however a few 
properties of the overall reactions involved: 

a. For describing the previously assigned changes in molar numbers 
BNj the transformation frequency 5Z j of the l-th overall reaction 

depends on the choice of the other overall reactions although the 
Individual overall reactions are thermodynamically independent of each 
other since they correspond to n independent equilibrium conditions. 

b. In particular cases from the velocities obtained for the overall 
reactions there may be derived the elementary reactions that are of most 
interest. 


■^See in particular the papers of the Bucken Institute especially 
in the Ztschr. physikal. Chem. B. since 193^* 

1 

Schafer, Kl., Ztsch. physikal Chem. B. B6, pp. 212-228, 19^0. 
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2. Relations Between the Overall and Elementary Reactions 

Among the three types of molecules Xq, Xg, and X^ which in our case 

represent three different excitation stages or tautomeric forms of the 
same chemical substance we have the following essential possibilities 


Reaction 1 
Reaction 2 
Reaction 3 



(105) 

(106) 
( 10 ?) 


They can be looked upon as elementary reactions or as overall reactions. 

In the latter case only two of these three reaction equations can be 
considered as thermodynamically independent. The frequency with which, 
for given molecular changes 8Nj, the Z-th reaction proceeds from left 

to right is as before assumed to be SZ In the case of equilibrium the 

velocities of the overall reactions will as_before be denoted by Uj, the 
velocities of the elementary reactions by u^ and u^, depending on 

) 

the direction. The questions as to whether in the equilibrium case u^ = uj, 

i.e. whether the principle of microscopic reversibility holds, may be left 
open. In Table 1 are given three possibilities for the choice of the 
overall reactions (denoted by the thick lines in the diagrams) and the 
relations between the 8N^ and the 8z 1 and between Uj and the 

various u^ and u m . For the 8z j and Uj the three possibilities 

of the choice of the overall reactions are distinguised by strokes above 
the symbols . 

Inspection of Table 1 leads to the following conclusions: Although 

In the two cases I and II the same overall reaction occurs, namely, 

X]_ { * X2, the frequencies 8zq and 8z]_ for describing any molar 

number changes 8Nj are in general different. For, on the one hand 
we have 8zq = 8Ng and on the other 5zi = -8Rq, and 8N2 and -BN-|_ 
need not be equal. Also the velocity frequencies Uq and Uq corre- 
sponding to the overall reaction in the two cases I and II are not 
necessarily identical. For it is not possible to derive from Table I 
that Uq = Up. There exist, however, between the various stroked Uj 
values certain relations which are given in the lower part of the table. 

They correspond except for the sign, which for the is always positive 

from the definition, completely to the relations between the various 
stroked Sz ^ . 
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Of the six different stroked Uj values shown in the right part 
of Table 1 only 3 are mutually independent and these can be obtained 
from sound dispersion measurements. For this purpose it is only necessary 
to compare the experimentally obtained dispersion curve with the two 
theoretically obtained curves for the overall reactions of case I and 
case II respectively. 

In the right part of Table 1 are found 6 equations with which for 
the equilibrium case the velocities u \ and tfj of the elemetary 

reactions are expressed in terms of the overall reaction velocities U m . 
Since altogether only 6 elementary reaction velocities exist (reactions 
to left and right counted separately) it might be thought that they 
could be computed individually from the U m values. In this way 
might perhaps be proved the principle of microscopic reversibility 
which in the case of equilibrium requires that the velocities of 
the elementary reactions to the right agree with the velocities 
of the elementary reactions to the left. This however is only a pure 
assumption which need not be satisfied for maintaining a state of 
equilibrium and which is obvious only for the over-all reactions in 
either direction. 19 


The six equations on the right of Table 1 for the and U i 

are not independent of one another. Writing them in the somewhat clearer 
form 


Ul 




u 2 

u 2 


— * 4— 

u 3 + u i 

— ^ 

u 2 + u 3 

f- f— 
U 2 + Ug 

— ) — ^ 

U 1 + u 2 


Ul 


ui + u 2 


(108) 


19t 0 assume the law of microscopic reversibility as a generally 
valid law of nature is equivalent to denying for the molecular aspect 
any type of check such as not only very frequently occurs in the 
microscopic aspect but is often found of very great use. I have in mind 
the "catch" type of control such as is used in rotating doors for 
regulating the traffic through them and which can also be observed in 
nature (for example for the artifically constructed cases of certain 
meat-eating plants). Naturally such a catch control in the molecular 
aspect can not be described by a fixed potential surface such as 
introduced for example by H. Eyring and M. Polanyi (Ztschr. physikal. chem. 
B. 12, 279 (1931) for the theoretical investigation of the activation 
heats. 
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The and uj are obtained as the quotient of two determinants which 

now both possess the value zero. The determinant in the denominator 
for example becomes 


1 0 
0 0 
0 1 
0 0 
1 1 



0 0 0 


1 



(109) 


It remains unchanged if to the first row are added other rows multiplied 
by the signs given on the right. There is then obtained a new zero- 
rowed determinant, so that the latter vanishes. For the corresponding 
numerator determinants the proof is quite analogous. 

If therefore, fo£ the equilibrium case it is desired to compute the 
velocities u^ and uj of the elementary reaction from the overall 

reaction velocities U^, assumed known, further assumptions must be made. 
The principle of microscopic reversibility might for example be 
tentatively assumed as correct and there would then be obtained 

u z = u z = u 2 ( 110 ) 


and therefore from equation ( 108 ) 


u x = iCU-L - U 2 + U X ) 

(111) 

U-2 = g-( - Up + U 2 + 

(112) 

u 3 = ^(u x + u 2 - Ui) 

(113) 


If in particular 


Ul = U 2 = U-l = U 


( 11*0 
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there would thus be obtained on the assumption of the validity of the 
principle of microscopic reversibility 


) 


u 


1 







(115) 


For the case of complete checking of the reaction in one direction, on 
the contrary, there would be obtained, as can be easily seen from the 
schematic diagrams on the left of Table 1, 


and 


— * (- X 

U 1 = u 3 = Ug = 0 


> 


i— 

u i = u^ = U2 


= U 


(116) 


or 



and 



(117) 


The question as to the manner in which the thermodynamic equilibrium 
between the three forms X-^, X 2 , and X^ is maintained, whether microscopically 

reversible or irreversible, remains entirely open and there exists as yet 
no experimental possibility of providing a reliable decision. Never- 
theless in very many cases the principle of microscopic reversibility will 
have to be assumed although on the basis of molecular theory the checking 
by means of the potential barriers in the sense of Eyring-Polanyi could 
be considered. Under these conditions it is possible from sound dispersion 
measurements for the equilibrium case to determine not only the 
velocities Uj of the overall reactions but also the velocities of uj 

of the elementary reactions. The given procedure, which was here 
carried out only for a special reaction kinetic system, can readily be 
generalized to other cases. 


3. High-Temperature Systems (Dissociation of 0 2 and C0 2 ) 

We now come to the main consideration itself of the present 
paper, namely the kinetic investigation of very rapid homogeneous gas 
reactions in the high-temperature range where dissociation effects 
assume prime importance. The application of the formulas will be 
clarified with the aid of two examples, namely, the dissociation 
of pure 0 2 and pure C0 2 . 



In the dissociation of pure 0r>, which ve shall consider at I = 2600°K and 1 atm, we are 
dealing with a single reaction which again we denote hy the index 1. We have 


0 


2 


2 0 


W-l = - 121.61 kcal 
1 

-S2- = 0.0243 Atrn2 

w 

RP ~ " 23»56 


(118) 

(119) 

( 120 ) 


( 121 ) 


Type of 

Particle j Index pj 

o 2 1 0.9760 


~ V il 
Nj J 1 

1.025 1 



20 9 . 340 to 9.022 


0 



0.0240 41.66 -2 

1.0000 v i » -1 



4.965 


9.235 to 8.925 


> 


It = 1.274 to 1.286 


( 122 ) 


h - KV i - [* 



I- 1.274 + 0.274 (23.56) 
- 1.286 +0.286 (23.56) 

t. 


5.182 

5.452 


(123a) 

(123b) 


20 Both values are taken from Johnston and Walker. The higher value from the Jour. Amer. 

Chem. Soc. 57 > 682 to 684 (1935)> the lower value from Jour. Amer. Chem. Soc. 55? 172, 187 (1933). 
The computation is here carried out with the two values, which differ appreciably, in order to 
learn the effect of such uncertainty on m and k - m. The higher Cp value is probably the 
more correct since for it the ^Ag level (7881.6 cm - *-) of the O2 is taken into account. 
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B 


11 




v x 2 = 23.56 + 1 (1.025) + 4 (hi. 66 ) 


1 = 190.23 
(124) 



Vl 


23.56 + 1 - - 22.56 


(125) 


Substituting these numerical values in equation (29) there is obtained 
for the differential isentropic exponents 


0 

Ai 


r 


A 1 

B 11 

- = < 

1.274 - 0.0874 = 1.1866 

(126a) 

1 

C 1 

Al 

B 11 

1.286 - 0.0949 = 1.1911 

(126b) 


Although in the example under consideration only about 1.2 percent of the 

0 2 is dissociated the difference between the isentropic exponent m and 
the ratio of the specific heats K is already quite considerable, namely 
O.091 + 0.004, i.e. 7.1 percent of the K value. 

We shall now estimate in what frequency range the sound dispersion 
by the 0^ dissociation is to be expected and what values of the sound 
velocity and sound absorption come into question. The excitation of the 
vibration heats will first be assumed to follow completely the sound 
frequencies, i.e. we shall for the present dispense with the Kneser case 
of the sound dispersion. 

In equation (70) we now require as a new constant cpqq, i.e. we 
must first estimate the reaction velocity Uq: 

3x 10 9 number of the bimolecular collisions per second of an 0 atom 
at 273° K and 1 atm. with other particles 

3 x lo9 (\/273/2600) = 1 X 109 a number of bimolecular collisions per . second 

of an 0 atom at 2600° K and 1 Atm. total pressure with other 
particles 

0.024 x 10 9 number of bimolecular collisions per second of an 0 atom 

with other 0 atoms at 2600° K and 1 Atm. total pressure. 
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0.024 x 10 ^ (lO ” 13 / 10 - 9 ) = 2.4 x 10 ^ SI number of collisions per second 
of an 0 atom at 2600° K and 1 atm. total pressure -with other 
0 atoms in the presence of a third particle of the collision 
complex of the two 0 atoms possesses a life of 10 “13 seconds. 

If each of the last-named termolecular collisions is suitable for 
a recombination of the 0 atoms there follows from equation (68): 

cp 11 = (fttf* 2 ) = 0.2l8f (127) 

2.4 X lo3 X 0.024 


The 2 in the numerator compensates exactly the error arising in the 
denominator from the fact that in the summation over all successful 
collisions each collision, on account of the two 0 atoms that take 
part, was counted twice. 

The order of magnitude of cp]_i can be estimated, however, not only 
from the recombination of the 0 atoms but also from the dissociation 
of the 0 2 molecule, though less accurately. Assuming a monomolecular 
reaction 

-122 x 10 3 

H 2 Z = 0.976 X 10 13 x e ( 128 ) 

N 


in which the factor lO^ -3 is uncertain by two powers of ten but still 
possesses the correct order of magnitude, then we have for 2600° K: 

— = 0.976 x 10 13 x 10" 10 - 25 = 5.49 x 10 2 (129) 

N 


and therefore 


<Pll 


2fff 

5.49 X 10 2 


0.0115f 


(130) 


The agreement between the two entirely different estimates which 
lead to (127) and (130), respectively, is quite satisfactory and can- 
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not be expected to be better . We continue the computation with the 
following numerical values: cp-Q = 0.2f and the following constants 
obtained from k = 1.27k: 


A x = 5.182 B ir = 190.23 C x = - 22.56 (131) 

Substituting in equation (70) there are obtained for the various 
frequencies f the numerical values given in Table 2 and figs. l(a) 
to 1 (c) . 

In figures 1(a) to 1(c) a logarithmic scale was chosen for the 
frequencies f along the axis of abscissas as is customary in repre- 
senting sound dispersion measurements. The magnitude m, the sound 
velocity a, and the amplitude damping 7 referred to unit length 
with increasing f pass through a point of inflection which lies at 
about 1^50 cycles per second. The limiting values of m, as was to 
be expected, are, for very low frequencies, the differential isentropic 
exponent for complete equilibrium attainment, and, for very high 
frequencies, the ratio of the specific heats k = 1.274-0 at constant 
pressure and constant volume. The sound velocity a changes in the 
entire dispersion region by about 3-6 percent. The amplitude damping 7 
per centimeter amounts to about 2 x 10"3 cm“^ at the point of inflection 
so that the sound intensity decreases over a distance of 10 centimeters by 

(l - e -27 x 10) ~ 27 x 10 = 0.04- times the initial intensity, or a lowering 

by h percent. 

The phase angle cp and the amplitude damping yl per wave length l 
with increasing frequency reach a maximum which similarly lies at 
about 14-50 cycles per second, the numerical value of cp being very small 
since cp = 0.0349 corresponds to 2°. The previously given approximation 
equations (78) to (80) are thus entirely 'justified and the term in 
parenthesis proportional to cp 2 can even be neglected. The maximal yl 
values here computed on the basis of dissociation are about 10 times 

as large as the former yl values which, for moist Og at room 

temperature, were found in about the same frequency range and are 
conditioned by incomplete vibration excitation. 2 1 

As a second high-temperature example we consider the dissociation 
of pure COg between T = 2000 and 3200° K at a total pressure p = 1 atm. 

Here we must consider a total of 4 types of particles, namely, 

CO2, CO, 2>2) and 0 which we denote by the particle index 1 to 4 and 
which are obtained on the basis of two overall reactions. The corre- 
sponding stoichiometric coefficients v are given in table 3- 

^Kneser, H. 0., Ztschr. techn. Physik 16, 2l6 (1935 ) > fig. 5* 
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There will again first he computed the isentropic exponent m for 
complete attainment of equilibrium. The required numerical values are 
given in Table 4. Since according to the latter the partial pressure 
of the 0 atoms is still quite small, it is natural to use two types of 
computations for the isentropic exponent m, namely, the exact computation 
taking into account the two reactions 1 and 2 (see Table 3 ) and an 
approximation computation in which only reaction 1 is taken into account 
and reaction 2 neglected. The isentropic exponent ' obtained by the second 
method will be denoted by m’. We thus have: 


0 

A 1 

A 2 

A 1 

B 11 

b 12 

Ag 

221 

b 22 

1 


V 

Cl 

B 11 

Big 

c 2 

^1 

b 22 


and 


0 

A 1 

A 1 

B 11 

1 

A 1 

0 

H 

B 11 


(132) 


(133) 


Substituting the numerical values of Table 4 there is obtained Table 5 
and figure 2. . 

The ratio K of the specific heats rises continuously in the 
temperature range 2000 to 3200° K because by the increasingly stronger 
dissociation the number of the 2- and 1-atom molecules increases at the 
cost of the 3-atom C0g molecules. K - m and k - m r are positive 
and increase with increasing temperature on account of the increasing 
dissociation. The isentropic exponents m and m' respectively run 
through a minimum because above 2600° K with increasing temperature k 
increases more rapidly than k - m and k - m' . The values m and m' 
obtained with and without account taken of the 0g dissociation 
respectively differ only very little. 
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The fact that the approximate expression (133) yields practically 
the same numerical value as the exact expression (132) is due to the 
circumstance that B22 is considerably greater than the absolute values 
of the remaining terms of the two determinants. If the latter are 
developed in minors there is obtained 


and 


A 1 B 11 
a 2 B21 

1 A-l 

Cl Bn 
C2 b 2 i 


Ag 


0 

Al 

b 12 

S B22 

A 1 

B 11 

B22 




a 2 


1 

Al 

B12 

a B22 

Cl 

B 11 

B22 





from which m JB m'. The expression 


(13*0 


B22 




2 N v 2 

5J 8 


is considerably greater than the other B^ because, on account of the 
stoichiometric coefficients v ji (see Table 3)? it contains a term N/N^ 

referring only to the 0 atoms and, since the latter is inversely 
proportional to the small partial pressure p Q of the 0 atoms, it is 
very large. 

We shall now consider the sound dispersion which can be expected 
in pure CO2 on the basis of the dissociation equations 1 and 2 
(see Table 3) at 2 600° K and 0.9994 atm total pressure. For this 
purpose we must first estimate the overall reaction velocities Up 
and Ug for the equilibrium case which according to equation (68J also 

determine the values cp^ and cp22 entering equation (70). 

We start with the reaction velocity U2 since it - is of small 
effect for the process of the sound dispersion and we can dispense with 
its consideration briefly. For if at complete attainment of equilibrium 
of reaction 1 the reaction 2 were no longer able to follow the increasing 
sound frequency f then according to Table 5 the m value of equation (JO) 



NACA TM 1268 


^3 


would in order of magnitude increase only from m = 1.1003 to m' = 1.1022, 
that is very slightly compared to the stage k - m' = 0.0728 which will he 
passed through if the reaction 1 does not occur together with 2. But even 
if the reactions 1 and 2 in the same frequency range can no longer follow 
the sound frequencies there is still no basis found for a larger effect 
of the reaction 2 on the dispersion process. We estimate as before the 
recombination of the 0 atoms in the termolecular collision: 

1 X 1 0 9 ~ number of the bimolecular collisions per second of an 0 atom 
with other particles at 2600° K and 1 atm total pressure. 

O.OO67? X 10^ = 6. 79 x 10^ X number of bimolecular collisions per second 
with other particles of an 0 atom in dissociating COg at 2660° K 
and 1 atm total pressure. 

6.79 X 10^ x 10 = 6.79 X 10 2 ~ number of collisions per second 
of an 0 atom in dissociating COq at 2600° K and 1 atm total 
pressure with other 0 atoms in the presence of a third particle. 

If each of these termolecular collisions leads to recombination 
there is obtained 


CP22 


2mf2 

6.79 X 10 2 x 0.00679 


2.72f 


(135) 


To estimate the reaction velocity U]_ we similarly use the reverse 
reaction, i.e. the combustion of CO with $2 at 2600° K in the presence of 
excess CO2 and the 0 atom content of the equilibrium. As elementary 
reactions that determine the velocity we shall take the following into 
considerations 22 

Elementary reaction I: 

CO + 0 + M > CO2 + M + ca. 125 kcal (136) 


Elementary reaction II: 

CO + 0g -■ ) C0g + 0 + ca. 8 kcal (137) 


For the reaction I we estimate: 

7 x 109jy number of bimolecular collisions per second of a CO molecule 

at 273° K and 1 atm. total pressure with other particles. 

— — - 

Jost, W. , Explosions- und Verbrennungsvorgange in Gasen, Berlin 

1939, P- 338. 
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7 X 10^ y273/2600 = 2.27 x 1C)9 ft* number of bimolecular collisions per 
second of a CO molecule at 2600° K and 1 atm. total pressure 
with other particles. 

2.27 x 10 9 x O.OO679 = 1.54 x 1 0 7 « number of bimolecular collisions 

per second of a CO molecule at 2600° K and 1 atm. total pressure 
with other 0 atoms. 

1.54 x loT x 10 -1 3 x 2.26 x lo9 = 3.5 x 103 2 number of collisions per 
second of a CO molecule at 2600° K and 1 atm. total pressure 
with 0 atoms in the presence of a third particle. 

If the elementary reaction I , which is assumed to occur for each 
termolecular collision, determines in general the velocity for the 
overall reaction 2C0 + O2 — > 2CO2 there would be obtained from 
equat ion _ ( 68 ) 


cp 


11 


2tt £ . x 5 — — ■ = 2.2 x io- 2 f 

3.5 X 103 X 0.1628 


(138) 


In the high-temperature range the elementary reaction II of 
equation (137) could also be an important factor. It undoubtedly 
possesses a considerable, though as yet unknown activation energy q-j- j . 

We shall estimate with various arbitrarily assumed qpp values and obtain: 

2.27 x lO^ x O.O78 = 1.77 X 10® ~ number of the. bimolecular collisions 

per second of a CO molecule at 2600° K and 1 atm total pressure 
with O2 molecules. 

The conversion probability per second of a single CO molecule 
with an O2 molecule under the same temperature and pressure conditions 
then becomes 


-20000 


1.77 X 

10 8 X e 

1.986 x 2600 _ 2.69 X 10 ^ 

for 

%I = 20 

kcal 



- AQQQ.Q 




1.77 X 

10 8 X e 

1.986 X 2600 _ 7.69 x 10 ^ 

for 

q n = kO 

kc al 



-60000 




1.77 X 

<D 

X 

CO 

3 

I.986 x 2600 _ x lo3 

for 

a =60 

H II 

kcal 



NACA TM 1268 


45 


If the elementary reaction II of equation (137) were alone to 
determine the velocity of the reaction 2C0 + Og — ^ 2CO2 there would 

he obtained 


9 


. £ jtf „ X_2 = 2.09 X 10" 5 f for q = 20 kcal ( 139) 


11 3.69 x 10 6 x 0.1628 


cp 

y ll 


= 1.005 X 10" 3 f for q_^ = 40 kcal (l40) 


7. 69 x 10 4 x 0.1628 


cp = 2£f_X_2 1 _ k.82 X 10 -2 f for a = 60 kcal (l4l) 

11 3 „ 11 

1.604 x 10 x 0.1628 


These values of cp computed on the basis of reaction II extend over 

the range 2 X lCT^f to 5 X 10 -2 f which includes the value 
cp n 1 = 2.2 x 10“ 2 f obtained from equation (138). This is 22 times as 
great as the value of cp i;L that is obtained from equation (l40) for 
the reaction II with 9^ = 40 kcal. As long, therefore, as q^ is 

not considerably greater tnan 40 kcal, it must be possible on the basis 
of sound dispersion measurements to distinguish between the two reaction 
possibilities I and II and at the same time it should be possible from 

tests at various temperatures to determine q_^. The relations are 

numerically even more favorable, i.e. q^ becomes also greater than 

40 kcal, if we drop two of the restricting assumptions used thus far, 
namely, that the reaction I is successful for each three-particle 


2 3case II, i.e. the reaction equation (137) is the more interesting 
as it takes place between two quite 'stable initial molecules and still 
leads to an active particle (0 atom) so that chains can be initiated. 

The knowledge of the reactions between such stable initial molecules 
where high activation energies must always be expected is quite gener- 
ally important for many starting processes (induction periods) such as 
for example occur in the prereactions in the knocking engines. 
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collision and that for the kinetic process of II only two quadratic 

ok 

energy terms are of' importance. If in the first case only the 
nth termolecular collision leads to a transformation and if, in the 
second case 6 quadratic energy terms were effective 2 ? there would be 
obtained in place of ( 138) : 


9-q = n X 2.2 X 10 -2 f (142) 

and in place of equation (l40): 


m = 1.005 X 10~ 3 f 

y ll 

1_ I 40,000 ' 

T3 \1.986 x 2600 

This means that the experimental differentiation between the reaction 
possibilities I and II would become still easier because the already 
large value (138) would become still larger and the small value (l4o) 
still smaller. 

We shall now see in what frequency range the sound dispersion in 
pure C0 2 at 2600 °K and 1 atm. total pressure is to be expected and 

what numerical values are assumed by the sound velocity and sound 
absorption. We shall compute the dispersion curves from equation (70). 


1.005 X 10~3f 

30 


(143) 


24 

The probability that a system with 
possesses a total energy greater than q 


s quadratic energy terms 
is given by 


5 


W s(q) 


Zl 

e* T X 


1 

r s /2 



zl 


so that for s = 2 
See H. J. Schumacher 


RT 


there follows the expression used V 2(q) ~ e ~ 

, Chemische Gasreaktionen, p. 15; Leipzig 1938. 


25 — 

4 terms from the vibrations of the CO and 0 g terms from the 

relative motion of the colliding molecules; or 2 terms from the 0g vibra- 
tion and 4 terms from the relative motion: see J. J. Schumacher loc. 


cit. p. 20. 
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The required values of Aj, sn<3. ^ l are gi ven in table 4. For 

<Pli are used the rounded numerical values of the expressions (135)> 
(138), and (l4o), that is. 


Case I: cpqq = .02f 

Case II: cp-j^ = .OOlf 

The results of the computation are given 


922 = 2f 
cp 22 = 2f 


(144) 


in table 6 and figures 3a to 3c. 


The m and a curves in figure 3a again show an S-shaped trend 
as was found to be the case in figure la for the O2 dissociation. In 
figure 3a for the case II two points of inflection can clearly be seen: 
the first at about 500 cycles per second for the small stage of the O2 
dissociation ( O2 2 0) and a second one at 1.5 X 105 cycles per second 
for the large stage of the CO2 dissociation 2C02<j^*2C0 + O2. 


The cp and yl curves possess in figure 3t marked maxima (compare 
also fig. lb for the dissociation of pure O2 at 2600° K and 1 atm.) 
which moreover in the two cases I and II are at the same height and 
possess approximately the yl value 0.1. The maximum yl values are 
thus evidently to a large extent independent of the assumed reaction 

velocity Uq for the CO2 dissociation.^^ 


The equally large maximum values of yl (damping coefficient per 
wave length) bring about the result that the damping coefficient y per cm 
in case II assumes on the average larger values than in case I. More- 
over in the reaction-kinetically determined sound dispersion only the 
yl values run through a maximum with increasing frequency while the 
7 values approach a constant maximum value, 27 which lies higher the 
higher the frequencies in the sound dispersion range) 7 max is ^ ca se II, 
for example, 20 times as large as in case I; see figure 30- 


^It may be noted that the maximum 
lies at about 0.1. 


7 l 


value In figure lb likewise 


^The reaction-kinetically determined sound absorption corresponds 
in optics to a- continuous absorption which sets in at a definite fre- 
quency range with constant intensity and which extends up to infinitely 
high frequencies. The case of light absorption in a finite frequency 
range restricted at each end, such as occurs, for example, in color 
filters, appears to have no counterpart in sound absorption". 
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In order to obtain a still clearer picture of the magnitude of the 
sound absorption determined by pure reaction kinetics there are shown 

in figure 3h the curves of e"-^ and e - ^?. These give the trans- 
mitted fraction of the sound energy if the sound wave! traverses a dis- 
tance of 5 and 10 cm respectively. In case II for 10? cps only 2.7 
or 0.074 percent of the initial sound intensity remains. Above 5. 10 5 cps 
no sound measurement is any longer possible if the sound must traverse 
a path from 5 to 10 cm from the sound source over a medium uniformly 
heated to 2600° K. It is not possible however to dispense with such 
measurement if sound of definite and controllable frequency can be pro- 
duced at 2600 o K. As a result of this too high sound absorption the 
experimentally accessible sound dispersion range is as yet in many cases 
for high-temperature reactions limited toward the high frequencies. We 
shall return to this question later. But first there will be discussed 
other effects on the sound dispersion which for reaction-kinetical investi- 
gations can similarly restrict the useful frequency range for the eval- 
uation of sound tests. 


Translated by S. Reiss 
National Advisory Committee 
for Aeronautics 




TABLE 1 


TRANSFORMATIONS AND EQUILIBRIUM REACTION VELOCITIES BETWEEN THREE DIFFERENT EXCITATION 
STAGES OR TAUTOMERIC FORMS OF THE SAME CHEMICAL SUBSTANCE 



Conclusions: Bzg = ES'j U 2 ** ^ 

8z^ = BZ^ % = U 2 
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TABLE 2 

ESTIMATED SOUND DISFERSION IE THE O 2 AT 2600 °K AND 1 ATM. 
• ON THE BASIS OF DISSOCIATION ALONE 


0-3 ' 

m 

q> 

a 

£cm s"^j 

_ 7 

[cm - 5 

7 l 

10 

1.1866 

4.80 x 10" 4 

9.006 X 10 h 

1.675 x 10‘ T 

1.508 X 10" 3 

1 x 10 2 

1.1869 

4.78 X 10" 3 

k 

9.007 x 10 

I.665 x 10" 5 

1.499 X 10" 2 

5 X 10 2 

1.1953 

2.155 X 10“ 2 

9.040 x 10^ 

3.746 x 10“^ 

6.773 X 10" 2 

1 X 10 3 

1.2132 

3.299 X 10 -2 

9.107 X 10^ 

1.138 x 10~ 3 

1.036 X 10” 1 

2 X 10 3 

1.2423 

3.400 X 10 ” 2 

9.217 x 10 ^ 

2.319 x 10 -3 

1.068 X 10” 1 

1 X 10 1 *' 

1.2721 

1.031 x 10~ 2 

9.325 x 10 4 

3.474 X 10“ 3 

3.239 X 10~ 2 

1 X 10 

1.2740 

1.054 X 10~ 3 

9.332 X 10 4 

3.550 x 10 -3 

3.312 X 10~ 3 
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Table 3 


REACTION EQUATIONS AND STOICHIOMETRIC COEFFICIENTS 
FOR THE DISSOCIATION OF C0 2 




TABLE 4 


vn 

ro 


NUMERICAL VALUES FOR THE COMPUTATION CF THE DIFFERENTIAL ISENTROPIC EXPQN 0 ST n 
IN DISSOCIATING C 0 g (l atm = 760 Eg) 


T to °K 

2000 

2300 

2600 

‘ 

2900 

3200 

22 ( 5 **.*.... 

PCOg 

0.002337 

0.01150 

0.0605 

0-2235 

o. 64 o 

— 1r> Atml /2 , . , , . 

0.000708 

0.00525 

0.0243 

0.0835 

0.220 

PCO2 111 Atm - 

0.9960 

.9100 

.7518 

.5185 

- 

.2912 

p^O In Atm ....... 

.01528 

.06045 

.1628 

.3106 

.4443 

in Atm 

.00760 

.03000 

.0780 

.1393 

.1760 

P 0 In Alan 

J 

p = y* ~ Pj Iji Atm ■ • • « 

.00006 

.00091 

.00679 

.03116 

.0923 

1 . 0189 ^ 

1.0014 

.9994 

.9 996 

1.0038 

K ........... 

1.1620 

1.1635 

1.1750 

1.1982 

1.2350 

Wp In fecal 

- 132.79 

-231.96 

-131.16 

-130.31 ■ 

-229.44 

VTp In fecal ....... 

-221.15 

- 121.41 

-121.61 

-221.78 

-221.90 

A 1 * * 

4.221 

3.558 

3.270 

3.286 

3.554 

*2 

3.749 

3.182 

2.945 

2.992 

3.275 

B u 

437.29 

131.96 

67.09 

49.40 

47.90 

B 22 

-101.62 

-5.50 

11.587 

14.453 

13.666 

B 21 * * 

- 204.55 

-7.80 

9-737 

12.973 

22,486 

B 22 * 

6.812 x 10^ 

4463 

624.16 

155-73 

67.39 

=1 

- 32.43 

-27.88 

- 24.40 

-21.63 

- 19.37 



- 29.50 

- 25.58 

- 22.55 

-20.15 

- 18.19 
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COMPARISON 


T In °K 

K 

K - m 
m 

K - m' 

m* 


TABLE 5 


OF THE DIFFERENTIAL ISENTROPIC EXPONENTS m AND m' WITH 
THE VALUE FOR DISSOCIATING C0 2 


2000 

2300 

2600 

2900 

3200 

I.l 6 l 0 

1.1635 

1.1750 

1.1982 

1.2350 

.0312 

• 0555 

.0747 

.0939 

.1195 

1.1298 

1.1080 

1.1003 

1. 1043 

1.1155 

.0310 

.0548 

.0728 

.O 896 

.1082 

1.1300 

1 . 1087 

1.1022 

1.1086 

1.1268 











» °K AHD 1 ATM. C®f 3 HB BASIS OF MBBOCIAIIOHS 
°2*=i2 0 


v n 


n 

1 

M 

e -10r 

e-207 


A.O67 x 10-3 

768.2 

1.0000 

1.0000 

2.906 X ID' 2 

76.95 

.9962 

. 992 * 

M? fll O u 1 A "2 
^OAU A XV 

-»n cm 

An/*« 

• yocK 

•9 7^6 

7.260 X U 3 - 2 

25.75 

.9721 

.9*50 

9*516 X 10 " 2 

15 . 5 * 

.9*07 

.88H9 

1.011 x io -1 

11.16 

. 913 * 

. 83*3 

9.558 X IO " 2 

7.853 

.885* 

.7839 

6,396 X ID -2 

3.956 

.8507 

.7237 

l .*35 X 10" 2 

• 7937 

.83*6 

.6966 

1.UA3 X 10-3 

.079* 

.8339 

.695* 

I.H3 X ITT 1 * 

.0079K 

-8339 

. 695 * 


1 . 5*7 x 10-3 

768.5 

1.0000 

1.0000 


2.81* x 10-3 

76.90 

.9996 

.9992 


t Vt/i^ *. 

X*M<J A iW - 

r-w / AA 

f *uyv 

•9315 

■9633 


*.07* X 10 -2 

2.567 

.8533 

.7281 


7.210 X 10* 2 

1.287 

.5710 

.3260 


9. *6 X 10 -2 

.777 

.2958 

.08750 


9.505 x 10-2 

-393 

.08901 

7.923 X 

10-3 

f, oUr v irb -2 

1 /VTfi 
**-7 1 ^ 

nV( 3 li)i 

■y u 

A. 

irt.? 
XV “■ 

2.008 X IQ -2 

.0793 

.0290* 

8.*33 x 

10 -* 

2.861 x 10-3 

.0079* 

.02719 

7.393 x 

10-^ 
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(a) Dependence of ra and. the sound velocity a on the fre-luency. 



(b) Dependence of the damping constant y referred to unit length (cm) for the 

sound amplitude. 



(c) Dependence of the phase angle <p and of the damping constant yl referred to 
the wave length l on the fre.juency. 


Figure 1.- Estimated sound dispersion in pure 0 2 at 2600° K and 1 atmosphere 
total pressure on the basis of the dissociation 0„^2 0 alone. 
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Figure 2.- Differential isentropic exponent in pure CQ? for full equilibrium: 
c p 

k = = ratio of molar heats, m = exponent in the isentropic curve 

pV m = const, taking account of the two equilibria 2C02* = ^2C0 + Og and 
0 2 ^ 2 0. m' = isentropic exponent for the case of the first reaction alone. 
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(a) Dependence of m and the sound velocity a 
the frequency. 


™/s 

790 

785 

780 

775 

770 

765 


on 



(b) Phase angle 9 81111 damping constant yl referred to 

the wave length Z as functions of the frequency. 


Figure 3.- Estimated sound dispersion in COg at 2600° K and 1 atmosphere 
total pressure on the basis of the dissociations 2CC>2 5^ 2C0 + 0g and 
0 £ <— 2 0 . 











